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Abstract 

^fs We present a scheme for full quantum state tomography tailored for two spin qubits in a double 

^^ quantum dot. A set of 15 quantum states allows to determine the density matrix in this two-qubit 

(N 

^ space by projective measurement. In this paper we choose a set gained from mutually unbiased 

.^ bases. We determine how those 15 projections can be represented by charge measurements after 

^^ a spin-to-charge conversion and the application of quantum gates. The quantum gates include 

^ exchange-based gates as well as rotations by electron spin resonance (ESR). We assume the ex- 

pi I perimental realization of ESR operations to be more difficult than the exchange gate operation. 

^ Therefore, it is an important result that the ESR gates are limited by a 7r/2 rotation for one of the 

^ electron spins per measurement. 
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I. INTRODUCTION 

Since electron spins in quantum dots were proposed^ as qubits, a lot of experimental 
progress has been made with a two-electron double quantum dot as the basic cell of this kind 
of quantum computing: Time evolution due to exchange interaction was demonstratecP as 
well as electron spin resonance (ESR) of one electron spinP Recently, Brunner and coworkers^ 
demonstrated both, exchange interaction and ESR rotation in the same double quantum 
dot. In principle this allows to perform arbitrary quantum gates in this system. Ultimate 
control of a quantum system can be demonstrated in experiment by state tomography, 
i.e., gaining full information of a quantum state, in general described by a density matrix. 
Since the density matrix is Hermitian with trace 1, it can can be described by A^^ — 1 
real parameters for a A^-dimensional Hilbert space. Therefore, from repeated projection 
of the unknown quantum state onto at least A^^ — 1 known states, the density matrix 
can be reconstructed from experimental data. A higher number of different measurements 
can increase the precision.l^ The experimental values of measurements are not perfect due 
to external noise and because they are not performed infinitely often. This may lead to 
unphysical results for the estimation of the density matrix by violating the non-negativity 
condition. By applying maximum likelihood™^ or Bayesian^HlS methods those results can 
be avoided. 

State tomography is well established in nuclear magnetic resonance experiments with 
qubits encoded by nuclear spins of trapped ions^^iiS for up to eight qubits.'^ It was also 
performed with superconducting qubits.'^^'^ For spin qubits, state tomography was done 
for qubits represented by the singlet and one triplet level of two electron spins in a dou- 
ble quantum dot for the single-qubilpSl and the two-qubit case^^l as well as recently for an 
exchange-only qubit built by spin states of three electrons in a triple quantum dot.l^ 

Here, we present a measurement scheme for determining the unknown mixed quantum 
state of a two-qubit system realized by two electron spins in one double quantum dot. A 
complete set of measurements for the two-qubit space is, e.g., given by James et al.,'^in that 
case referred to the polarization of two photons. We use another series of projections, which 
is based on a set of so-called mutually unbiased bases^SHlll and which includes entangled 
states in contrast to Ref. [3 In this issue our scheme also differs from the scheme used in 
Ref . 1211 where the state of each of the two qubits, both encoded in one two-electron double 



quantum dot, is measured separately and the density matrix is determined using correlation 
data. The measurements we propose are projections on one-dimensional subspaces of the 
four-dimensional Hilbert space, and thus joint measurements on the full two-qubit space. 
An important result of our paper is that these measurements can be done by applying ESR 
maximally for a 7r/2 rotation on one of the electron spins per measurement. We assume the 
ESR to be experimentally more costly than exchange gates as in experiments the time for 
a complete ESR rotation was in the order of 100 ns^ and a complete rotation due to the 
exchange interaction could be achieved within some ns.'^ Furthermore, we assume the ESR 
gates to have a lower fidelity than exchange gates. 

In Sec. [Tljwe describe the double-dot system with two electrons and the possible manip- 



ulation of the quantum states. Sec. Ill provides information on the matrix space for the 
traceless part of two-qubit density matrices and on properties of a basis in this space rep- 
resented by pure quantum states. The main result of this paper, a scheme for constructing 
such a set of quantum states using quantum gates and spin-to-charge conversion, is given in 



Sec. IV In Sec. ^ we discuss how the influence of imperfection in the quantum gates might 



be treated in experiment, before we conclude in Sec. VI 



II. PHYSICAL SYSTEM 

We consider a double quantum dot with two electrons taking only the lowest (orbital) 
energy level in each dot into account. Therefore, due to the antisymmetry of the wave 
function (Pauli principle), we have a six- dimensional Hilbert space since we have one state 
with both electrons in the left dot, one state with two electrons in the right dot, and a four- 
dimensional space with one electron per dot. We will denote charge states with rii electron 
in the left and n^ electrons in the right dot with {ni,n^). The Hamiltonian, taking into 
account the difference between the energy levels, e, the Coulomb penalty for both electrons 
in the same dot, U , and the hopping matrix element between the dots, t, reads 

H = i^{ni-n2) + -^ fiiifii - 1) + t ^ (c^ci, + c\/2s) (1) 

j=l,2 s=t,4 

where the operators of the occupation number in the dots, -hi and n2 are given by hi = 



^s=t, 



cLci. 



;H 



We will consider the states with charge configuration (1,1) as logical Hilbert space and 



regard the electron spin in each dot as a quantum bit. The aim of this paper is to give a 
scheme for determining an unknown density matrix of this two-qubit space. In the hmit 
|t| ^ |t/ ± £:| the (1,1) states are approximately decoupled from the states with charge 
configuration (2,0) and (0,2). For this case the effective Hamiltonian for (1,1) states can be 
obtained by applying a Schrieffer- Wolff transformation,'^^! 

H.. = J "^"7"' = -JPs (2) 

with J = At^U /{U'^ — e^), and where ctj is the vector of Pauli matrices for the electron spin 
in dot i = 1,2. Note that H^s is given by —J times the projection operator Ps on the spin 
singlet IS") = "^ because only an antisymmetric state can couple to the (0,2) and to 
the (2,0) space. The strength of this exchange coupling, J, can be tuned in experiment by 
varying e or t. 

We now consider in addition to the exchange interaction above a magnetic field yielding 
the Zeeman Hamiltonian 

i^B = hi • CTi + h2 ■ cr2. (3) 

Let us restrict the discussion first to a magnetic field only in 2;-direction, hi = /ii^e^, 
h2 = h2zGz- The homogeneous part of the magnetic field, h^ = ^^i+^^2 ^ gpjjts off the 
triplet states, because |T+) = |'|''|') and |T_) = \ID obtain the Zeeman energy ±2/i^ whereas 
l^o) = _^ is not affected. A gradient in the magnetic field, described by Ah^ = hiz — h2z, 
leads to a mixing of IS*) and \Tq). If \Ahz\ and \hz\ are large compared to At'^/U, which is 
the value of J for e = 0, then the energy eigenstates with charge configuration (1,1) are 
approximately |T+), |T_), \'\i), \it)- At higher values of e, \S) becomes an energy eigenstate, 
see Fig. [l] Different directions in the magnetic field in the left and the right quantum dot 
can additionally induce an avoided crossing between |'|''|') and \S). 

Projective measurements on the two-qubit space are done by charge measurements after 
a spin-to-charge conversionj ^^ * ^^ ' which we describe in the following. For e = the four 
lowest energy levels belong to states with charge configuration (1,1), which constitute the 
logical subspace of the system. The starting point is an arbitrary state \ip) in this logical 
subspace and e = 0. By changing gate voltages, e can be increased to a value higher than 
U, modifying the spectrum to a situation where the state (0, 2) is the ground state of the 
system, see Fig. [I] (a). Performing this varying of e backwards in time would map the (0,2) 
charge state on one specific state in the logical subspace, lips)- The forward procedure thus 
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FIG. 1: (Color online) (a): Spectrum of a two-electron double quantum dot with a magnetic field 
with \hz\ ^ IA/12I ^ At'^/U, hz, A/i^ < for small values of e and e close to C/. The avoided 
crossing between IS") and |tt) ^^ the blue dotted line is induced by slightly different directions of 
the magnetic fields in the left and the right dot. (b): Three ways to tune e to achieve projection 
on Itl), dash-dotted (green) line, on \S), dashed (red) line, and on |tt), solid (cyan) line. Note 
that the e-axis and the time-axis are broken 



converts the {ips) amount of {ip), given by {i'sli') li^s), to the (0,2) charge state and the 
remaining contributions of \tp) are still in charge configuration (1,1). Therefore, measuring 
a charge state (0,2), corresponds to a projection onto the state {ips) in the logical subspace. 
Now, we will consider three different cases shown in Fig. [I] (b) and explain which is the state 



\ips) for the different time dependences of e. Assuming we are in the situation shown in Fig. [I] 
(a), then a sfow change of e from to a value higher than U would be an adiabatic transition 
of the |T+) state to the (0,2) state. If we increase e slowly from up to a value below the 
avoided crossing between IS*) and |T_|.) and then proceed fast through this avoided crossing, 
the charge readout effectively projects onto the state |ti)- An overall fast transition yiels a 
projection onto \S). We will use this three different kinds of projection in our scheme for 
state tomography in Sec. IV In the state-tomography experiment for the 5* — Tq qubits^^Ell 



the projections onto |ti) and onto \S) were used. To allow projections onto other than the 
three states described so far, we assume quantum gates to be applied prior to the sweep. In 
principle more complicated transitions are possible, they are described by a Landau-Zener 
transition.'^^'^ State preparation is done in the opposite direction by converting the (0,2) 
state in a well-defined (1,1) charge state. 

It has been shown that the exchange interaction together with arbitrary single-qubit 
rotations allows for universal quantum computing,'!' i.e., every unitary operation can be 
constructed from these elementary gates. To achieve full control on a single qubit, a tunable 
magnetic field in one direction is necessary with individual tuning for each dot. Additionally, 
rotations on the Bloch sphere around another axis than the 2;-axis must be possible, which 
can be done using ESR. This requires an oscillating magnetic field orthogonal to the axis 
of the time-independent magnetic field. The frequency of the oscillation is given by the 
Rabi frequency of the constant field. As a sufficiently strong high-frequency magnetic field 
is hard to realize experimentally, Brunner et al.^ used an inhomogeneous magnetic field and 
achieved an effective time-dependence in the Zeeman Hamiltonian by applying an ac voltage 
shifting the electron back and forth. Due to an offset in Zeeman energy between dots it was 
possible to perform ESR on each of the dots separately. Since exchange interaction was 
performed as well, in this double dot sample, operations allowing for full control on the 
two-qubit space have been demonstrated experimentally. 

III. GENERAL CONSIDERATIONS ABOUT TWO-QUBIT TOMOGRAPHY 

We want to find a scheme for state tomography for the two-qubit space. This means the 
scheme has to provide the full information of the density matrix p, which is a Hermitian 
4x4 matrix with trace 1. Therefore, we need at least 15 different measurements to deter- 



mine the 15 real parameters describing p. In the space of 4 x 4 matrices we can introduce 
the scalar product {A\B) j^,^ := ti{A'^ B) and the basis {0^='^-^-'^^ with A;=4i + j and i,j G 
{0, x=l, y=2, 2;=3}} where o"ij := ctj ® 1 and 0"2j := 1 (8) cr^ are Pauli matrices for the first or 
second spin, and cr„o = 1- Note that this basis is orthonormal with respect to (•|-)a,/ because 
{aii(72j\o'ikO'2i) j^j = 'iSij6ki. We can expand p in the {Dk} basis, 

15 

p = Y,PkDk- (4) 

fc=0 

We know already that all pk are real because p is Hermitian and that po = 1/2 due to the trace 
condition. The aim of the tomography is to determine every parameter p^, fc = 1, . . . , 15. 

Estimations for pi, . . . ,pi5 have to be acquired from experimental data. In the physical 
setup we consider in this paper, measurements are done by measuring the charge state after 
increasing the detuning e. This converts a state from the logical spin space (qubit space) with 
charge configuration (1,1) to a (0,2) charge state. In the previous section we found conditions 
under which this leads to a projective measurement on either \tD, \S), or |r+). We assume 
that universal quantum gates are feasible and therefore, any other state could have been 
mapped by a unitary operation onto the state which is converted to the (0,2) charge state. 
Therefore, spin-to-charge conversion and universal quantum gates allow for a projection 
onto an arbitrary state 10). Each measurement within the tomography scheme will project 
onto a state \(f)j) and thus yields tr(|0j) {4>j\p)- The projection operator Pj = |0j) {(j)j\ is 
also a Hermitian 4x4 matrix with trace 1, which means {Pj\Dq) ^^ = ^ (Pj|l)^^ = 1/2. 
Nevertheless, Pj can be expressed in terms of the six real parameters^ for the state |0j). 
For Pj given in the ^^J^^ basis, see Appendix [Aj We are seeking a set {|0j)} that provides 
full information of p. If 15 projectors fulfill the condition that {Pj — 1/4, j=l, . . . , 15} is 
a basis in the space of traceless 4x4 matrices, we call these projectors a minimal set or 
quorum. 

We introduce the invertible 15 x 15 matrix V with 

V,k = {Pj\Dk)M with J, /c = 1, . . . , 15, (5) 

i.e., considering only the traceless contributions of the projectors Pj. Since tr (Pjp) = \ + 
^k=i'PjkPk, the parameters pk are given by pk = ^j=iV^j^{tT{Pjp) - |). We obtain an 
experimental estimation rrij for ii{Pjp) by dividing the number of the charge state (0,2) 
outcomes by the number of experimental runs for the jth projection experiment. This 



statistical data does not yield the exact values of pk but estimations for it, which we call pk 
and which are by linear reconstruction given by 
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Although the maximum likelihood method would be preferred over the linear reconstruction 
for practical purposes we want to consider Eq. ^ in order to obtain information about 
the influence of V on the statistical error expressed by the covariance matrix Cki = E{{pi — 
Pi){pk— Pk)) where E{-) is the expectation value of the probability distribution, which is given 
by the density matrix p and the set of measurements Pi, . . . , P15. Note that E{-) is in general 
not expressed by the quantum mechanical expectation value tr (Op) of an operator O, but 
this is the case for those operators which are directly measured, i.e., for every j = 1, . . . , 15 
we have a binomial distribution with the number of trials, Nj, and the probability of success 
tr(Pjp). The probability distribution for the indirectly obtained stochastic variables pk has 
to be gained using the matrix V. Using Eq. ([6]), we find for the covariance matrix 

(15 15 \ 

j=l i=l / 

(7) 

= E ^fc"/ ^(K- - tr(P,p))K - tr(P.p))) V^\ 



i,j=i 



-Bji 



which is equivalent to C = V^^B(V^^)^. Our aim is to find a set of measurement for which 
every value \Cki\ is small. The covariance matrix of the measurement results, B, is given 
by Bji = (5jj tr(Pjp)(l — tT{Pjp))/Nj < 5ji/{ANj) as a property of the binomial distribution 
describing rrtj, and the upper limit is reached if tr(Pjp) = 1/2. Using the adjugate matrix, 
{B:d]V)ki = (— 1)*^^' det(P///,fc) where V/i/k is the matrix one gains by deleting the /th row 
and the fcth column of P, to express the inverse of V we obtain 



IG 



kl\ 



f. 1 itr(P,p)(l-tr(P,p)) 



J = l 



A^,- 



15 



|(adjP)fc,(adjP); 
4A^jdet(P)2 



M 



(8) 



A 



The length in matrix space of P,- — 1/4, i.e. a row of V, is \ {Pj — \\Pj ~ |) 
determinant is invariant under the Gram-Schmidt orthogonalization, yielding upper limits 
for |det(P)| as well as for the matrix elements of the adjugate matrix of P, which are 
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determinants of 14 x 14 matrices gained from V by deleting one row and one column. For 
adj V we find |(adj V)kj\ < (f )^- Assuming Ni = N2 = ■ ■ ■ = Ni^ = N, this results in the 
following inequality 

\r I < ^^(i)'' ~ 0-06682 

' '^^ - 4iVdet(P)2 ^ Ndet{Vy ^ ' 

We now concentrate on maximizing the value det(P)^. The upper bound for |det('P)| is 
given by (|)^ ~ 0.1156, which would be reached if the rows of V were orthogonal, which 
means the 15 projectors Pj have to fulfill 

(P, - l/4|Pfc - 1/4)^, = ^ {Pj\Pk)M = I {<Pj\<Pk) r = 1/4 for J ^ k. (10) 

In Appendix [B] we prove that this is not possible. 

The quorum from Ref. [7] for the two-qubit space, after translation from photon polar- 
ization states to spin states, reads {|tt), \U), lit), lis,t), \iyi), 114), Itxt), \txiy), \txtx), 
lU.), in.), i;t.), in.), in.), lU.)} where It,) = ^, \iy) = Mi), and |t.) = ^1 . 
We have removed the state \ID because 15 states are sufficient here. These basis states are 
all separable and lead to | det('P)| = 5J2 ~ 2 ■ 10"^. 

A better result can be obtained by the well-known concept of mutually unbiased bases.'^SESl 
Two bases in Hilbert space are called mutually unbiased if a state of the first basis has the 
same overlap with all states of the second basis. The idea of mutually unbiased bases 
originally referred to a situation where A^ -|- 1 observables are measured,'^ and each of this 
measurements yields the probability for being in one of the corresponding A^ eigenstates, 
which are here chosen to be basis states of one of the mutually unbiased bases. In that 
case it was shown that the mutually unbiased bases represent an optimal set to determine 
p in the sense of minimizing the statistical error.EZI For the projective measurement by 
the spin-to-charge conversion, which we consider in this paper, the situation is slightly 
different since every projection is done by another measurement. In particular, it is not 
necessary to measure all 20 probabilities; it suffices to perform the projective measurement 
for three states out of the four states in each of the five bases, providing a basis in the 
15 dimensional matrix space, and consequently a minimal set of projectors. It has been 
shown that for a A^-level system a set of A^ -|- 1 mutually unbiased bases exist if A^ is 
the integer power of a prime number.'^ Hence, for our four-dimensional case, five bases 

{|0O,l) , 100,2) , 100,3) , 100,4)}, • • • , {|04,l) , • • • , |04,4)} fulfilling 

I (0jfc|0/m) P = 5j74mH y^- (H) 



can be found. Note that states from different bases with numbers j 7^ / fulfill | {(pjk\(pim) 



|2 



|, which is the condition (10). This means that the matrix space spanned by the traceless 
parts of the corresponding projectors contains five three-dimensional orthogonal subspaces. 
On the other hand, the states in the same basis are orthogonal in Hilbert space, {4>jk\4'jm) = 
5mk- We show in Appendix [C| that this always leads to | det(P)| = ^ = 0.03125. 

IV. EXPLICIT CONSTRUCTION OF SPIN-QUBIT QUORUM 

We use a set of mutually unbiased bases given in Ref. |2H] and present an explicit con- 
struction of the following 15 states out of these bases, forming a quorum, 

1 + 0"l2 + (72z + (7lz(72z 



^1) = Itt) , Pi 

^2) = in) , A 



4 

1 + 0"l2 - (T2z - (^lz<^2z 



4 

^3) = SWAP Iti) = lit) , Ps = ^^^ + "" " ""^^' 



4 

_ 1 + (^Ix + 0"2x + Crix'72a: 

4 — 



V'4) = e*^'^^^v/SWAPe^f'^-|5) = |t.tx), A ^ 

^5) = e*2-i^ 1^4) = |;^|^) , P5 = , 



^7) = e-f (--+-^) 1^4) = \tyty) , i^7 



4 

4 ' 



^s) = e^^^^MV^r) = lU.) , P3 = i^i^l^.±^^.^f^l^, (12) 

V^9) = e^^-^ iV^r) = lU.) , A = ''^"^^"7""^^"^^ 



^12) = e'^'^'^ iV'io) , P] 



12 



4 
1 + ai^a2x + cria;a2y - 0-1^0-2^ 



4 



^14) = e^^-^ iV'is) , Pi4 = i^^l^^ - ^^^^^^ - ^^-^^^ 



^15) = e^^-^^^ iV'is) , Pi 



15 



4 

1 - aiya2x + 0"l2Cr2y + (7lx(72z 



4 

where the states \ip3i+i) , |'?/'3J+2) , iV^Sj+s) are belonging to the same basis, i = 0,1, 2, 3, 4. The 
first three bases {i = 0, 1, 2) are given by the product of spin states in each dot where the 
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FIG. 2: The quantum circuits used to realize the quorum (12). Note that ESR is used in the 



preparation of the states IV'4) and jV'io), the states IV's) , • • • , {tpg) and jV'ii) , ■ ■ ■ , jV'is) can then be 
gained by rotations due to a static magnetic field from 1^4) or iV'io) respectively. The gate -Rxf 
on the upper line in panel (b) denotes a rotation of the first qubit around the x-axis by the angle 
7r/2, i.e., the operation e*4°"i^. The other rotations Rx and Rz are defined accordingly. The angles 
61 and 02 in panel (d) depend on the index j of the final state but they are the same for the maps 



1-04) I— 5- IV'j) and iV'io) I— ^ IV'j+e)- The actual values of 9i and 02 can be found in Eq. (12), e.g., for 
j = 9 they are given by 9i = —tt/2 and 02 = vr/2. 

quantization axis is chosen to be z, x, and y, thus those states are separable. The other two 
bases (z = 3, 4) consist of maximally entangled states. 

Within the scheme we are using projective measurements on the states |tt)) Iti)? and 
IS*) as described in Sec. [Ill In order to access the states in the set above not directly given 
by one of these states, we apply various quantum gates directly available for spin qubits, 
see Fig. [2] The SWAP gate is defined by SWAP |siS2) = \s2S1), i.e., it interchanges the 
information of the first and the second qubit. It can be realized by applying the exchange 
interaction (2) for a time te given by ip = Jq" J{T)dT = it, where we neglected the magnetic 
field (gradient) which is possible if J is large compared to the gradient Ahz- Consequently, 
the a/swap gate can be achieved by applying the exchange interaction for yielding (^ = 7r/2. 
Since a^ + 022 commutes with the exchange coupling, the infiuence of h^ can be treated 
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formally as if it was applied before or after the exchange interaction. Furthermore, we need 
an ESR pulse on the first spin to access the states \ip4) , . . . , iV'is)- We want to emphasize 
that for each of these twelve states we need only one ESR 7r/2 rotation, e^i,"^^ of only one 
of the spins. As we assume the fidelity of ESR operations to be lower than that of exchange 
gates or rotations around the 2;-axis, this is a crucial point for an experimental realization 
of a tomography scheme in a double quantum dot. 

Note that ESR cannot be completely avoided because exchange interaction and rotations 
due to a static magnetic field together with the projections on W\), \S), |tt) can only 
access the five-dimensional subspace in matrix space spanned by {<yizi<y2zi<^iz<^2zi<^ix<^2x + 
criycr2y, (yixO'2y — criyCr2x}- This cau be seen by regarding the time evolution under exchange 
interaction, which leaves Pft ~ -^i ^^^1 Ps invariant and evolves P^i = P2 to 

i^Ps p -i^Ps _ 1 + (^iz(^2z + cost^(ai^ - a2z) + smcp{a^^a2y - aiy(T2x) ,-. ^. 

and the time evolution under a gradient Zeeman field in z-direction, leaving Pf^ and P^ 
invariant and evolves Ps to 

e r-se ^ 

(14) 
The single-qubit gates containing aiz or 0-22 can in principle be achieved just by waiting 
for a certain time with negligible J and an applied magnetic field. Since the magnetic field 
cannot easily be switched off completely this requires precise control of the time which is 
spent between state preparation and the projective measurement. Considering for example 
the states 1^4) and {ipr), we see that the state with spins polarized in x-direction, {ip^), and 
the one with both spins in y-direction, {ipr), are connected just by a rotation around the 
2;-axis for both spins, which is realized by a (large) constant magnetic field. But clearly it is 
necessary to distinguish those states \1lj4) and \ipi) experimentally to get the full information 
of the density matrix p. Therefore, state tomography requires more control than would be 
necessary to demonstrate ESR rotations. Note that the adjoint quantum gates included in 



(12) are applied on the initial state in the tomography experiment. Afterward a lack of 
time control might not harm the experiment as a e*'^'^°"i^"'"°'2z) g^^^g results in no more than a 
global phase on the states \ti), Itt)? and \S). Only the projection on these states by the 
spin-to-charge conversion has to be precise. 
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V. ESTIMATING THE ERROR 

Measuring a quantum system always produces statistical results, thus the exact quantum 
state cannot be determined by an experiment which is performed a finite number of times. 
But increasing this number of repetitions cannot lead to an arbitrary high accuracy of the 
estimation for the quantum state either since the measurement scheme unavoidably contains 
systematic errors. Therefore, it is necessary to know how precise the desired measurement 
scheme is actually performed. We assume the error in the measurement setup to arise 
mainly from the applied quantum gates, and among those the ESR gate, 6*4^^^=", to have the 
largest infiuence an the total error. In the experiments described in Ref. H] the fidelity of 
the exchange rotation could not be extracted because two ESR rotations were also included 
in the explored two-qubit gate affecting strongly the fidelity of the combined gate. Errors 
in the quantum gates may occur due to fiuctuations in the magnetic field, for instance, 
caused by nuclear spins,'^'^!! -vv^hich is most damaging for the ESR gate. Charge fiuctuations 
mainly infiuence the exchange interaction and thus the SWAP and a/SWAP gate. Finally, 
differences in the times for which interactions are actually switched on, affect all quantum 
gates but in particular the spin rotations around the z-axis according to a static magnetic 
field. 

Knowledge about the experimentally realized precision of the projective measurements 



suggested in Sec. IV, can practically only gained by prior experiments. One option to 



achieve this is a self-consistent tomography^SES] -^j^grg ^i^e perfect states from Sec. IV , P. 



3 



(j = 1, . . . , 15), are replaced by imperfect and unknown ones, P/. This implies that before 
using the set of measurements to determine an unknown quantum state, the set {-Pj} itself 
has to be determined first. For this purpose, one could, e.g., perform measurements to 
obtain statistical knowledge on the quantities 

M,, =tr(/^'P;), ^,j = l,...15. (15) 

Note that preparing the system in the state P[ is the inverted procedure of the projective 
measurement on this state, staring from a (0,2) charge state. The states P' can be determined 
approximately afterward via a maximum likelihood algorithmP^'Sa ^^p ^q a basis in Hilbert 
space which has to be fixed as Mij is invariant under a unitary transformation. Note that 
here, less information is needed than for quantum process tomography.'^ 
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Since the 15 states in Sec. IV are generated by a smaller number of different quantum 
gates it might be an alternative to describe those gates with a set of parameters including 
mean fluctuations and determine these parameters experimentally. In general, the imperfect 
projection operators Pj would be assumed to depend on a set of n parameters cti, . . . , «„, for 
a parametrization of fluctuating magnetic and electric fields see Ref. 1251, for a parameter fit 
of a ESR experiment see Ref. [31 For experiments which are repeated often, we can consider 
the mixed states P' calculated from the statistical distribution of the parameters ai, . . . , a^ 



) "n; 



Pj = dai...danPi{ai)-..Pn{<yn)Pj{<yi,...,an). (16) 

For Gaussian distributions it is sufficient to determine the mean values of a, and af, which 
means that n^ numbers have to be extracted from experiments. Note that errors in those 
values will directly lead to errors in the tomography of an unknown state, but it is not 
necessary that the distributions are sharp. Nevertheless, the deviations of the mixed states 
P' from the pure states Pj worsen the result of the tomography measurement as illustrated 
by the following example. Assume the matrices P4 — 1/4 and Pg — 1/4 to have the same 
directions in matrix space as P4 — 1/4 and Pq — 1/4, i.e., they are given by 



1 - /2 4/2 - 1 

3 3 



P; = -—^1 + A^i^,, (17) 



where fj = tr ( \ufPjP'-^,fPj\ is the fidelity of the quantum state Pj with respect to Pj, 
and where we assume /4 = /e = / for the reminder of the discussion. For a completely 
mixed state P', = 1/4, we find / = 1/2. Whereas p^ from Eq. dll) in the ideal case is given 

2(tr(P4p)+tr(P6p))-l 
P4 = ^ , (18) 

for the mixed states it is represented by 

3[2(tr(Pip)+tr(P^p))-l] 
P^ 2(4/2 _ 1) • (19) 

For the statistical error of p^ desired to be smaller than 6, the statistical error of the ex- 
perimental values for tr(P4p) and tr(Pgp), which we assume to be equal, has to be below 
S' = {.^ ■ A Chernoff bouncP^J^Sl p^n be applied to estimate the number of experimental 
runs needed to realize a desired precision.'^SESl xhe possible outcomes of one experimental 
run are for no projection and 1 if a (0,2) charge state is measured, the expectation value 
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is tr (Pjp). Repeating this experiment A^run times thus leads to a binomial distribution. The 
probability for the experimental value, i.e., the number of runs with result 1 divided by A^runj 
to be out of an interval [tr(Pjp) — (5',tr(Pjp) + 6'], Pout, is limited by 

Pout<2exp -^ . (20) 

V -' 'run / 

Therefore, the number of experimental runs which provide a Pout below a desired limit Pi 

can be chosen by 

^ ln(2/P,) _ 91n(2/P0 

This means that the reduced fidehty in the projection states can be compensated by increas- 
ing the number of experimental repetitions. 

VI. CONCLUSION 

We have presented a measurement scheme by accessing the states of a quorum for state 
tomography in the two-qubit space by realistic quantum gates and three types of spin-to- 
charge conversion in a two-electron double-quantum dot. The quantum gates include SWAP 
and vSWAP gates realized by exchange interaction, as well as ESR rotations. We assume 
the ESR gates to have the largest impact on the fidelity of the quantum states. Therefore, 
it is important that the ESR rotation is applied for no more than a 7r/2 rotation per state. 
Moreover, it should be an advantage that only one of the spins has to be treated with ESR 
since in the experiment by Brunner and coworkers,^ the fidelity of ESR was different for the 
left and the right quantum dot. 

Whether state tomography can be realized experimentally in the system which we consider 
here depends on the realization of the different quantum states. A reduced fidelity in the 
experimental representation of the (pure) quantum states does not prevent the tomographic 
reconstruction in principle but reduces its quality. The crucial point is the knowledge of the 
experimenter on which (mixed) states the unknown quantum state is actually projected. A 
lack of information about this state will in any case limit the precision of the tomography. 
Self-consistent measurements^^'^ or a calibration by analyzing the tomography dataP^ can 
reduce this lack of knowledge. In general, the accuracy of a realized tomography experiment 
has to be determined experimentally as well. Our measurement scheme is only one possible 
solution for state tomography in the two-spin system. The question which scheme is optimal 
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regarding the accuracy of the estimated density matrix with or without a hmited number 
of experimental runs remains open. The solution will crucially depend on the fidelity of the 
different quantum gates performed in the system. 
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Appendix A: Representation of a projection onto a state in matrix space 

Assuming that we can perform any unitary operation in our four-level system, we can 
project on any state 

10,) = «, in) + bj Iti) + 9 lit) + d, III) (Al) 

where |ajp + |6jp + |cjp + |(ijp = 1, without loss of generality, Qj G IR>o as a global phase 
does not matter. 

The corresponding projector can be written in the matrix form for the basis { |tt)) \ti)j 
lit), lit) }, 



/ |ajp Ujb* ttjC* ajd* \ 



ha* 16 P he* hd* 

* /-,* I 1 2 J* 

Cja^ CjOj \Cj\ CjUj 

\ djtt* djh* djC* \dj\^ J 



i,l=0,x,y,z 



(^li(^2l 



(A2) 



with 



16 



n 



jOO 



l«.f + l&.f + |c,f + M.f _i 



rijOy 

fT-JTO 



n 



■jxy 



Hi. 



jyo 



n 



wy 



njzo 



n 



■jzy 



2 

a*bj - b*aj + c*dj 


2' 
- d*c, 


2i 
a*Cj + b*dj + c*aj 


+ d*bj 


2 

a*dj - d*aj - b*Cj 


+ c*b, 


2i 
a*Cj + b*dj - c*aj 


- d*h 


2i 
b*Cj + c*bj - a*dj 


- d*aj 


a,p + |6,2-c,p 


-d,^ 


2 


+ d*c, 


2i 





a*bj + b*aj + c*dj + d*Cj 
njox ~ 



Hi, 



jOz 



'''jxx 



n 



]XZ 



f^jyx 



rijyz 



"^jzx 



"^jzz 



|a,f 


- \b, ' 


2 


- \d, ' 


a*dj 


+ d*a^ 


2 

+ b*c, 


+ c*b, 


a*c, 


- b}d, 


2 


- d]h 


a*dj 


- d*a, 


2 

+ ^.*9 


- ^'h 


a*c, 


- m 


2i 


+ ^i^J 


a*b, 


+ b]a, 


2i 


- dY^ 


l«jf 


- \h ' 


2 

- r- 2 


+ M.P 



(A3) 



Appendix B: No orthogonal basis in traceless matrix space by projectors on single 
pure states 



For the reconstruction of the density matrix from the measurement data, it would be 
ideal to have a set {Pi, . . . , P15} which fulfills for i ^ j 



P, 



^ I M 



1 1 



(Bl) 



4 '■ '^-^'-"' 4' 

which means that the traceless parts of the projection operators Pj would be orthogonal 
with respect to (■|-)a^. Analogously, for a single spin, the density matrix can be determined 
by projections on the x-, y-, and 2;-axis of the corresponding Bloch sphere. However, for 



our four- level system, Eq. (Bl) cannot be fulfilled for 15 projection operators of the form 



P,- 



) (0j I which we show in this appendix. For the proof we use the following lemma. 



"J/ \vj 



Lemma: Let V\ and Vi be two real Euclidean spaces with dimension n^ and n2. Let 
{ai, . . . , a„^+„2} be an orthogonal basis in the space Vi© V2 with |ajp = c > 0, \Py^aj'\^ = ci, 
and \Pv2Cij\'^ = C2 Vj = 1, . . . ,ni + 77.2 where Py^ (Pva) is the projection on the space Vi (V2). 
Then ci and C2 fulfill 

(B2) 



Ci _ ni 
C2 ^2' 



Proof of the Lem,m,a: We denote ai, . . . , a„^+„2 in an orthonormal basis 



{ei 



, e„^, e„j+i, . . . , eni+n2} where the vectors ei, . . . , e„^ are completely in Vi 
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and Cni+i, ■ ■ ■ ,eni+n2 are completely in V2. Now we can write the normalized vec- 
tors aiZ-y/c, . . . , a„^+„2/"\/c as the columns of a matrix A with the matrix elements 
^ij = (ei|aj)/v^. The matrix A is an orthogonal matrix, i.e., it fulfills A'^A = 1 but 
also AA"^ = 1, which refiects the fact that rows and columns denote orthonormal bases in 
Vi © V2. From the conditions in the Lemma we find 

m ni+n2 

^4 = !^ aud Y. 4 = 7. (B3) 

1=1 j=ni+l 

and thus 

ni+n2 ni ni+na ni+n2 

E E4 = (^1 + ^2)^ and 5^ 5^ 4 = (ni + n2)^. (B4) 

j;'=l i=l j=l i=ni+l 

On the other hand, using that the rows of A denote orthonormal vectors in ]R"'i+"2 yields 

ni ni+n2 ni ni+n2 ni+n2 ni+n2 

EE4 = Ei = ^^ ^^d ^ E4= E 1 = ^2. (bs) 

i=l j=l 4=1 i=ni+l j=l i=ni+l 

Building the ratio directly leads to 

j=l Z^i=l ^ij JB4} Ci |B5| ril _ .g„s 

ES"Er=:::i4 " ^^ " n,- ^ ^ 

Now we come back to the 15 basis states for the space of traceless Hermitian matrices 
constructed by pure states in four-dimensional Hilbert space. We use the notation from 
Appendix lAJ Assume that |(/)i) = |'|"t), if it is not, we can perform a unitary transformation 
to match this starting point. The corresponding projector is Pi = (1 + cti^ + 0"2z + criz<72z)/4:. 
It follows that aj = 1/2 in order to fulfill | {4>i\4>j) P = 1/4 for j = 2, . . . , 15. Therefore we 
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get 



njoo - 2 
Im bj + c*dj — d*Cj 



n 



n 



n 



n 



fT'jzO 



n 



JUJ/ 


ReCj 


2i 
+ b*dj 


+ d*bj 


jxU 


Imdj 


2 


+ <^]h 


j'-^y 


Imcj 


2i 
+ b]d, 


- d]h 


jyu 


2i 
b*c, + c*b, - 


- Re dj 


'jyy 
1 

_ 4 


+ \b, ' 


2 
- l9 ' 


-d,^ 




Imbj 


2 
- ^}d, 


+ d]c, 



jzy 



2i 



n 



jOx 



njoz 



n 



J XX 



n 



■jxz 



n 



■jyx 



n 



jyz 



n 



JZX 



n 



Rebj + c*dj 


+ d*Cj 


i- b, 


2 




Re dj - 


2 

f 6*c, 


+ ^b, 


Re Cj - 


2 

- b]d, 


5 

- d*b, 


Iradj - 


2 

f b*c, 


J 


Ira Cj - 


2i 
- b]d, 


5 

+ d*bj 


Re bj - 


2i 
- c*dj 


1 


\-b, 


2 

12 _ 1 


1 
-■ 2+ d-P 



■jzz 



Now we define a new basis of tlie space of 4 x 4 matrices, {tj; i = 0, . . . , 15} with 



T% 



ru 





1 




(^Ix + (^lx<^2z 


i'2 


2^2 




(^ly + 0"lj;CT22 


'4 


2v^ 




O'lxO'22; — ^l?/^2j/ 




2V2 




f^lx - 0-lxO"22 


'8 


2v^ 




0"!?/ ~ <^ly<^2z 


'iU 


2V2 




0"lxO'2x + 0"l?/0'2j/ 




2v^ 




0^1^ - 0-22 

^14= „ /K ' 



0-22 + O-I2 + 0'l20"22 

n = ^^;f , 



Ts = 

n = 

Tj = 

Tg = 
ri5 



2^3 

0'2x + 0'lzO'2x 

2^2 ' 

0'2y + <^lz<^2y 

271 ' 

0"la;0"2y + <^lyO'2x 

271 ' 

0'2a; — Criza2x 

271 ' 

0'2?/ — <^lzO'2y 

'~ 271 ' 

CI'lxCI'2y — (^lyO'2x 
~' 272 ' 

_ q'l2 + (^2z - 2<7lz<y2z 

2^6 



(B7) 



(B8) 



which again fulfills {Tk\Ti)j^^ = 6^- P2 — 1/4, ...,Pi5 — 1/4 are now supposed to be an 
orthogonal basis in the space spanned by r2, . . . , Tis. We show that this is not possible for 
the following reason. When we expand Pj in the new basis, Pj = '^i'lTijiTi, we find easily 
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for j > 2 from Eq. (B7) that 



^jo = 2' "^Ji = 0' 



m.2 



ReCi 

~7i 



-, "^i3 



mj5 = ^7=-, mje 



and thus 



E 

i=2 



"^,,, 



72 



2 _ l^jf + |Ci 



Re6j 



m 



, ,,Dj 



■j4 



m,7 



ImCj 

Imdj 

^2" 






1 — la,- 



On the other hand due to tr(Pj) = 1 and P^ = Pj, we find for j = 2, 



15 



15 



15 



15 



1 = {Pj\Pj)m = Yl ^% = ^% + "^|i +J2^% = j + Yl ^% 



i=0 



"^l 4 

i=2 i=2 



which yields 



15 

E 

i=2 



2 3 



Using Eq. (BIO) we find 



15 



i=8 



3 

4 



5Z'^?* = I-5Z'^^ 



i=2 



and thus 



E7 2 

i=2 ^ji 



M 



6 



(B9) 



(BIO) 



(Bll) 



(B12) 



(B13) 



(B14) 



V^15 2 

which means that the ratio of the projections of Pj on the subspace span{r2, . . . jTy} 
and spanjrg, . . . , Tis} is a fixed value which does not correspond to the dimensions of 
these subspaces and thus violates the lemma above. Therefore an orthogonal basis for 
span{ri, . . . , ris} cannot be constructed by Pi — 1/4, . . . , P15 — 1/4. ■ 



Appendix C: Reconstruction with mutually unbiased bases 



T^ if the matrix V as introduced in Sec. 



Ill 



Here we show that |det('P)| 
is gained from five mutually unbiased bases, including three states from each basis, 
{|0oi) ) 1002) , 1003) ) |0ii) ; • • • ; 1043)}- The rows of V represents the traceless parts of the 
corresponding projection operators P-^j+k = |0jfc) (0jfc|5 which we denote 



vi = m\D,) 



M 1 



{Pi\Di,)u)- 



(CI) 
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We want to apply Gram-Schmidt orthogonalization on the row vectors Vi which leaves the 
determinant of V invariant. As 

'Psj+k'Psi+i = I {(pjk\(pii) 1^ ~ 4 = % ( '^'^' ^ 4 ) ^^^^ ^5 J e {1, ... , 5} and k, I e {1, 2, 3} 

(C2) 
we only have to orthogonalize the five three-dimensional subspaces for i = j. From 
Vsj+k'Pjj+i = ^ki — \ "we gain easily the orthogonahzed row vectors from the Gram-Schmidt 
method, 

"P Sj+l = 'Psj+i, 

'P'3j+2 = Vs,+2 + lv3,+i, (C3) 

1/1 
'P'sj+l = 'Psj+S + ^V 3j+2 + g'Psj+l, 

which have the lengths -i/|, -i/|, and -i /I. A unitary transformation does not change the 
determinant and transforms V' to a diagonal matrix. Therefore, we find 

det(P) = det(P') = u/f • yi • yi) = (i) ' = ^. (C4) 
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